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The theory of thermodynamic properties of two-band superconductor with reduced 
density charge carriers is developed on the base of phonon superconducting mechanism 
with strong electron-phonon interaction. 

This theory is adapted to describe the behavior of critical temperature T^, energy 
gaps Ai, A2, and the relative jump of electron specific heat {Cs — Cm)/Cm in the point 
T = Tc along with the variation of charge carrier density in the compound MgB2 when 
substitutional impurities with different valence are introduced into the system. It is 
shown, that according to the filling mechanism of energy bands which overlap on Fermi 
surface, the quantities Tc, Ai, A2 decrease when this compound is doped with electrons 
and remain constant or weakly change when the system is doped with holes. The theory 
qualitatively agrees with the experimental data. 

Also is shown that the consideration of inter- and intraband scattering of electrons on 
impurity potential improves this agreement. 

1. Introduction 

The superconductivity in MgB2 compound with transition temperature ~ 39K 
was discovered about five years ago [1]. This discovery has generated a big interest and 
has attracted the attention of many researchers. Such interest is caused by big values of 
Tc, high values of upper critical field and critical current jc in this compound that put 
this relatively simply substance along with complex high temperature oxide compounds. 
As a results of many researches is clear that the superconducting properties of MgB2 is 
impossible to describe using the BCS - Eliashberg theory. The properties of this substance 
are characterized by a whole series of the anomalies which cannot be framed into this 
theory (see, for example, [2], [3]). The anisotropy of the system plays an important role 
in this case, particularly the overlapping of energy bands on Fermi surface which leads 
to the appearance of two and more gaps in the energy spectrum. The band structure of 
the compound MgB2 confirms such overlapping of energy bands [4], [5]. The possibility 
to describe this compound by applying the two- band theory of superconductivity was 
expressed by the authors of a number of works (see, for example, [6]). However, the model 



of superconductor with the overlapped on the Fermi surface energy bands was proposed 
long time before the discovery of high-temperature superconductors [7], [8]. The history of 
the development of the theory of multiband superconductors can be traced, in particular, 
in the surveys [9] - [12] . The Moskalenko's model for the two- band system [7] assumes the 
formation of Cooper pairs in each band and the tunneling of these pairs as a whole from one 
band into another. This model makes possible to obtain high values of Tc, two energy gaps 
Ai and A2 (in this case the conditions 2Ai/Tc > 3, 5 and < 3, 5 can be satisfied (in 

the one-band case 2A/Tc = 3, 5), the low values of the relative jump of electronic specific 
heat in the point T = (for MgB2 this jump is equal to 0.8, in one - band case - this 
is an universal value, equal to = 1-43), the anomalous temperature dependency of heat 
capacity in the superconducting phase, the positive curvature of temperature dependency 
of upper critical field He.-, and other properties [9]. The compound MgB2 possesses such 
properties, and the two-band model [7] describes the qualitative picture of the behavior 
of diverse physical characteristics of this compound. 

Together with the pure MgB2 the influence of chemical substitution of the atoms of 
boron and magnesium by other elements (for example, C, Lj, C„ etc.) on thermodynamic 
properties of MgB2 is of interest and is sufficiently intensively investigated [13], [14]. 
The thermodynamic properties of two-band superconductor with the lowered (variable) 
charge carrier density in the assumption of weak interaction in the case of non-phonon 
superconducting mechanism are investigated in the work [15]. A more general case, when 
all possible electron couplings (intra- and inter-band) are considered in the work [16], and 
also in the review [12]. 

At the present many works are published where the properties of the compound MgB2 
are investigated on the basis of the Eliashberg equation considering the presence of two 
energy gaps. Thus, this compound is examined from the point of view of strong electron- 
phonon interaction (see, for example, the works [17], [18] and references in them). 

The purpose of this work is to develop the theory of thermodynamic properties of 
superconducting states of the compound MgB2 with variable charge carrier density. The 
two-band model [15] and electron-phonon mechanism of superconductivity, inherent to 
this compound, are assumed. In this case the relative alignment of energy bands is 
considered, and the theory parameters are evaluated considering strong electron-phonon 
and Coulomb interactions. This approach makes possible to trace the energy bands filling 
and determine the influence on superconducting state in MgB2 of chemical substitution 
of magnesium and boron by other elements. At the same time the interband scattering 
of charge carriers on impurity potential is considered [19]. 

The work is structured as follows. In the Section 2 on the basis of Prohch's Hamiltonian 
the equations for temperature Green functions, diagonal M{x, x') and non - diagonal 
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E(a;, x') mass operators are obtained. The transition to band representation (nkQ) is 
performed expanding the Green functions and mass operators on Bloch functions. The 
system of equations for order parameters Ai and A2 with renormahzed constants of 
electron-phonon interaction Xnm (n;m = 1, 2) considering the dependency of mass operator 
Mn(k,n) from these constants is obtained. This system of equations is completed with 
an expression, determining the chemical potential /i. In the Section 3 the limiting cases 
T Tc and T = are examined. The basic equations of superconducting theory are 
defined in each of these temperature ranges. The analytical expressions for the critical 
temperature Tc, energy gaps Ai and A2 are obtained for specific range of values of charge 
carriers (chemical potential /x). The expressions for the jump of electronic specific heat 
Cs — Cn also are obtained in the point T — Tc. In the Section 4 the theory parameters 
for MgB2 compound are evaluated and the dependency of thermodynamic quantities on 
charge carrier density is built on he basis of energy bands filling mechanism. 

In the last section an additional mechanism of impurity infiuence on superconduc- 
tivity, namely interband scattering of electrons on impurity potential, also is considered 
[19]. The dependency of Tc value from the concentration of the introduced substitutional 
impurity C, Li and Cu is studied. The obtained results are compared with experimental 
data. 

2. The Hamiltonian of the system and the main equation 



The Hamiltonian of the system has the form: 

H^Ho + Hi (1) 

where Ho is the Hamiltonian of non-interacting electrons and phonons, and Hi - the 
Hamiltonian of electron-phonon interaction, determined by the expression 

Hi = gJ2[ dxij+{x)M^M^)- (2) 

Here g - the electron - phonon interaction constant, '4'^{x) and ipai^) - the operators of 
appearance and annihilation of electrons with the spin a in the point x, (fi{x) - phonon 
operator. 

Having using the diagrammatic method of the perturbation theory [20] for normal 
G{xx') and abnormal F{xx') temperature Green functions, we obtain the system of equa- 
tions 

G{x,x') ^G^{x,x') + j j dxidx2G'^{x,xi)M{xi,X2)G{x2,x')- 
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-J J dxidx2G^{x,xi)E{xi,X2)F{x2,x'), (3) 
F{xx') — j dxidx2Go{xi,x)'E{xi,X2)G{x2,x')+ 

+ J dxidx2G^{xi,x)M{x2,Xi)F{x2,x'), (4) 

where the sign ~ indicates the complex conjugation of the corresponding value, Gq{x, xi) 
- the electronic Green function for free electrons. 

The diagonal M(a;i,a:2) and non-diagonal T,{xi,X2) mass operators take the form 

M{xi, X2) = -g^D{xi, X2)G{xi,X2), (5) 

E(xi, X2) = -g^D{xi, X2)F{xi, X2). (6) 

Here D{xi,X2) - phonon propagator. 

Let's pass in the expressions (3) - (6) to the nkD. - representation, expanding the 
quantities from these formulas on Bloch functions [21]. Then, the system of equations for 
Green functions Gnmik, k' ,Q) and Fnm{k, k' and the expression for mass operators 
Mjnm'iP,!^, f2) and 'Lrnm'ip-iP' i ^) IS not difficult to obtain. Such approach makes possible 
to investigate the thermodynamic properties of two-band systems taking into account all 
possible intra- and inter-band electrons couplings. The limit of weak electron-phonon 
interaction in these systems was examined in the works [12], [16] . 

Here we will examine a much simpler case, considering only the approximation of 
Green's functions diagonal on bands indices. This case leads to the two-band model, 
proposed in the works [7], [8]: 

G'nin2(^l) ^2, = Gn-i^{k\^Vt\)5nin2^ki,k2^ 

{ki,k2,^l) = Fn^{ki,Qi)6 

nin2 

Mmm'{P,P,^) = Mm{p,^)Sj^jf6mm', 

^mm'{P,^, ^) = ^m{P, ^)Sp-p'Smm'- (7) 

In this approximation the solutions of the system of equations (3) and (4) for Green 
functions have the form 

G„(*,a) = -lM|I±£., (8) 

^"A^^-TtIS' ^"*^) = T#7!!' 
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where 

= 1 - ^ImM^ik, n), inik) = en{k) + i?eM„(fc, Q). (10) 
The expressions for mass operators Mn{p, fl) and E„(p, fl) can be written in the form: 

M^ip, f^) = E E E aLM h)D{p- h, n - fii)G„,(fci, fii), (11) 

E„(p, = E ki)D{p- n - nr)F^,{k„ n,), (12) 

where 

x{mpF, npp) ^ j df U^^^ (r) C/„p- (r) , (13) 

Unp(j) - the amphtude of Bloch function, O - Matsubara frequency for the electrons, Vq - 
the volume of unit cell. 

Let's introduce the definition S„(r2) = Z„(r2)A„(r2) and simplify the equations (11) 
and (12) disregarding the delay effects. After that let's integrate over Qi and then let's 
pass to the integration over energy according to the dispersion law of electrons energy 
from n - th band: 

2m„ 

Then, on the basis of relations (9)- (12) we obtain the system of equations 

Y th:'iEj2^ - th/3E2/2 ^ 

Ai = An J dei^^Ai + Xi2 J &2 ' A2, 

^ y, th(3EJ2 ^ - y, th(3E2/2 ^ 
A2 = A2i / dsi ; A1 + A22 / &2 ^J' A2, (15) 



2Ei J 2E2 

di —d2 



where 



-^nm = 9nm ^m, K = \J (^n " /^)^ + A^ , = — , A, 



A 



nm 
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dn^H-Cn, 4n = Ccn-/^, (n;m = l,2), (16) 

Njn - the density of electron states on m - th cavity of Fermi surface, Cm - the cut-off 
energy in the n-th band. According to the definition (10) and the dispersion law (14), for 
the quantities Z„ we obtain: 



Zi — 1 + All + Al2, Z2 — 1 + A22 + A; 



21- 



(17) 



Since we examine a system with the variable (including low) charge carrier density, 
let's supplement the system of equations (15) with a relation, determining the chemical 
potential /i (charge carrier density n): 



n 



E^{k) - \€m{k) - IJ,\ , 2\e„Xk) - fJ^ 



+ 



Em{k) l + exp(3Em{k) 



(18) 



When we have low values of charge carrier density n(// ~ A„) is necessary to correlate 
the solution of system (15) with the expression (18). 

Let's examine the phonon mechanism of superconductivity, which is inherent to MgB2. 
The cut-off of integrals in the integration over energy in each band are carried out on the 
characteristic phonon frequency = — Cn): 



<^0n 



(/in > UJon) 



dm. 



<^0n 
Cm A* 



(Cc 
(Cc 



1^ < l^On) 



(19) 



3. The limiting cases T ~ and T = 



Let's study the system of equations (15) at T ~ Tc. For this purpose we expand the 
parameters in terms of the difference (/9 — Pc), where P = 1/T: 



Cn(/3-/3c)^/^ + cW(/9-/9c)^/' + - 



3/2 



(20) 



Introducing this expansion into equation (15) and equalizing the expressions with the 
same degrees of the difference {P — Pc) , we obtain a system of equations which determines 
the coefficients and c^^ (n = 1, 2): 



Cl = -Aii/i(/i)Ci - Ai2/2(A«)c2, C2 = -A2li"l(//)Ci - X22h{lj)c2 



(21) 



Cl 



Pc 



,(1) 
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A 

~Pc 



e2ifi) - P'tF2il^)ci C2-Xl2Wc^\ 



^2 — 
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A 



Cl + 



/3c 



e2{ll) - l3tF{ll)ci C2 - A22/2(/^)4 



(1) 



(22) 



where 



•(") = -/ 



th(5ce/2 

2£ 



denize 



dx 



-dnl3c 



shx — X 



th — h th—— 



T ^cn T dji 



(23) 



Equalizing to zero the determinant D of the system (21), we obtain an equation for 
superconducting transition temperature Tc, and on the basis of (18) - an expression, 
which ties n and /i: 



n 



1 + Xuhin) + \22HlA + dh{ii)l2{lj) = 0, 

Yl [Ccn - Cn - \Ccn - fi\ + \Cn - ft 



where 



A11A22 ~ Ai9A' 



12^21- 



(24) 



(25) 



On the basis of the system of equations (21) is not difficult to determine the ratio of 
the parameters = A1/A2 at T — Tc. 



Ci _ Ai2/2(/i) 
C2 l + Aii/i(/x) 



1 + A22/2(/^) 
A2l/l(Ai) 



(26) 



The quantity Zc implicitly depends on charge carrier density (chemical potential) 
through the dependency of quantity Tc on this parameter, and also explicitly - due to 
the integration limits in formulas (23). 

The solutions of the system of equations (22) can be represented in the form: — 
Di/D, = D2/D where D, Di and D2 are the respective determinants of this system. 
The condition D — corresponds to the solvabihty of the system of equations (21) for 
the superconducting transition temperature. Consequently, it is necessary to put the 
condition Di — D2 — 0. From this condition the following relation results: 



Ni9M + N292{tJ,)/z\ 



PI N^F,{^^) + N2F2i^^)/zt^ 



(27) 
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For certainty let's assume Ci < C2 < Cci < Cc2- As it follows from the dispersion law 
of electrons energy (14), these parameters determine the mutual arrangement of energy 
bands. 

In the range of values dn/Tc and dcn/Tc ^ 1 the equation (23) allows to obtain 
analytical expression for superconducting transition temperature Tc- Considering cuqi — 
<^02 = and cutting off the integrals according to (19), on the base of (24) we obtain 

= ^e-^ (28) 

TT 

where 

_ All + A22 -a-f± y(Xn + A22 - o^y- - 4^ 

2a • ^^^^ 

Here the quantities (/? and ip are determined by chemical potential 

1. At we have ji' — (2 < (^o and — > (^o 

(p=-ln , ^ = l-A22-ln ; (30) 

Z (jJq Z Uq 

2. /i — > (^0 and Qi — P > ^0 

99 = 0, ^ = 1; (31) 

3. fl' — (2 > <^o and — < 

(^^lln^£l ^, -0 = 1 - Aii^ln^^i ^, fi^-^: j2'^jiZi/Z2. (32) 

Z (X^Q z (x^o -"^l 

In the points fJ. — C2 and = Cci we have respectively. 



2Aii + A22 ± Y (2Aii — A22)^ + 8A12A21 
= , (33) 



All + 2A22 ± y/(Aii — 2A22)- + 8A12A21 
" 2^ ^'^^ 

At Ci < < C2 and < ^ < Cc2 the value of Tc is determined by electron - phonon 
interaction in the first and second bands, respectively In this case the dependency from 
H takes place at p, — (i < coo and Cc2 — fi' < (^o- 



re = -V(/^-CiVoe-^/^?i, rc = -V(Cc2-/i'Vo e-^/^-. (35) 

TT TT 



Hereinafter AJi = An at A^2 = 0, and A22 = A22 at A^i = 0. 

For the jump of specific heat in the point T — TcWe have [15] 



Cs — Cn 



(3!^N^F^{l,)cl (36) 



V 

Using the relations (27) and (28) we bring the jump of specific heat to the form: 

Cs-Cn _ [N,eM + N^QM/^l? .^7^ 

In this case the relative jump of specific heat takes the form: 

Cs-Cm [iViei(/i) + iV202(/i)A 



If 



where 



Cn 4[iViFi(/.) + N2F2{ii)/zt][N^^i{ii) + iV2<^2(/^)] 

x^dx 



(38) 



J (l + e^)(l + e ^) 

-(M-Cn) 

At T = on the basis of given above expressions (15) and (18) is not difficult to 
obtain the system of equations, which determine the order parameters A^(0) = and 
chemical potential /x(0) = /i. This system takes the form 



n 



A. = VA^nA„/n ^":+V^- + ^- , (40) 
2V -dr. + Jdl + A2 



(41) 



E^n ken -Cn- ^ {Con - U? + A^ + V^(C, - /i)^ + 



In the analytical studies we will examine the specific ranges of values of chemical 
potential, when the following inequalities are fulfilled: A^/Jcn, A„/(i„ <^ 1. We use the 
definition of quantities d„ and den in accordance with cutting of integrals over energy at 
electron- phonon interaction (19). On the basis of the system of equations (40) in the 
interval where energy bands overlap (C2 < A* < Cci) is easy to get: 

a A2 
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Here Dq and z are determined by the value of chemical potential In the interval 
C2 < A* < Cci we obtain the following relations: 1) Jx' — C,2 < ojq (ci ~ ft > ^0, Dq — ujq, 

aln2; = All - A22 - Miz + — - In —\ ; (43) 

z 2 \ ujQ ) 

2) /i' - C2 > t^O, Cci- fl> i^O: Do = UJo, 

— — — 12 

a In 2; = All - A22 - A212; H ; (44) 

z 

3) /i' - C2 > '^o, Cci - < ^0, ^0 = \/^o(Cci -/i) 

aln;2 = All - A22 + — - A21-2 - a^ln^-^^^^5 — ■ (45) 

z 2 Cci - 

In the interval where the overlapping of energy bands is absent ^ < (2 and /x > Cci one 
energy gap is present Ai and A2, respectively. In this case for low (// — Ci < "^o) and high 
(Cc2 ~ A*') < values of chemical potential we have 

Ai = 2^cuo(/^-Ci)e-'/'", A2 = 2^a;o(Cc,-/^')e-'/'^^- (46) 

The formulae (46) show that the parameter Ai increases together with the chemical 
potential fi at the beginning of filling area of first energy band, and the parameter A2 
diminishes in the end of second energy band filling area. In the energy bands overlap- 
ping area C2 < A'- < Ca the picture is more complex and can be obtained as a result of 
numerical evaluation of given above expressions. In the case of low values of parameters 
dnidn ~ ^n) IS uecessary to obtain numerically a self- consistent solution for the system 
of equations (40) and (41). Let's note, that on the basis of these equations the analytical 
expressions for A„ and /i were obtained in the work [22] for the state of the deep Bose- 
condensation of local pairs (/x„ < 0, (A„//i„)^ <^ 1). 

4. Application of the model to the study of doped MgB2 superconducting 
properties 

The obtained above equations for T^. (24), (28) and for A„ (40) contain variable charge 
carrier density n (chemical potential /i) and, consequently, can describe the behavior 
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of the corresponding values of two-band superconductor when its atoms are chemically 

substituted by other elements. 

In particular, the compounds MgB2-xCx, Mgi^^ CuxB2 and other [13], [14] are of 
great interest. Numerous theoretical and experimental studies prove, that the anomalies 
of the observed physical quantities in the superconductible compound MgB2 can be un- 
derstood basing on energy bands overlapping on the Fermi surface and the presence of 
two energy gaps (see for example, [3] - [6]). Diverse approaches in scientific literature 
there are to study the compound MgB2. These approaches lead to ambiguous values of 
theory parameters. 

According to the contemporary researches [2] - [6], [18], [23], the compound MgB2 is a 
two-band superconductor with strong electron-phonon interaction. On the Fermi surface 
two energy bands overlap: a two-dimensional band S and a three - dimensional one vr. 
Hereinafter we will call the a - band the first band, and the vr - zone - the second band. 
In accordance with this let's designate the electron-phonon interaction constants An and 
A22 intraband, A12 and A21 - as interband ones. 

The table 1 gives the values of these parameters, and also the Coulomb interaction 
parameters /i*, obtained by the authors of works [18], [23] (a), [24] (b) and [6] (c) on the 
basis of Eliashberg equations. As it can be seen from this table, the values of the quantities 
Xnm and renormalized quantities Xnm aren't single- valued ones. A common thing for them 
is the relation An > A22, A12, A21 which tells about the significant role of two-dimensional 
a band in Cooper pair's formation process in MgB2. 

Solving the Eliashberg equations for the pure M gB2 on the basis of these parameters 
is possible to obtain values for superconducting transition temperature Tc, close to 40 K 



Here we put the problem to investigate the influence of substitutional impurity on 
thermodynamic properties of two-band superconductor. Since the impurity introduction 
cannot change the picture of phonon spectrum so, so that this change would radically 
affect superconductivity, to solve this problem we will use the given above renormalized 
equations of the two-band theory (24) and (40), in which the renormalized constants of 
electron-phonon interaction are determined by the relations: 



[18]. 



All 



A. 




A 




1 + All + A 



'22 



1 + A22 + A; 



12 — 



1 + All + A12 



12 



21 




(47) 



1 + A22 + A21 

where /x*^ is the parameter describing the Coulomb interaction of electrons. 
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The quantities \nm we determine from the condition that the values of Tc and Ai , A2 
would correspond to experimental data. 

We select the values = 39, AK, Ai = 7,lmeV, A2 = 2,7meV, 2 = |i = 2, 63, 
N1/N2 = 0, 73 and ljq = 75meV, which were obtained in the work [23]. The obtained in 
such a way the effective constants of electron-electron interaction are given in the table 1 
(case d). We use these values of to calculate basing on (31) the ratio Zc — A1/A2 and 
the relative jump of electronic specific heat (C5 — Ci\f)/Ci\[ (38) at T = Tc. We obtain: 
Zc = 3,2; {Cs — Cn)/Cn = 0,79. This, z < z^., and the jump of specific heat is small 
in comparison with the one-band superconductors (1-43) and is close to the experimental 
values [2]. 

Table 1 





All 


A22 


A12 


A21 






K2 


/^21 


All 


A22 


A12 


A21 


a 


1,017 


0,448 


0,213 


0,155 


0,21 


0,172 


0,095 


0,069 


0,362 


0,172 


0,054 


0,053 


b 


0,96 


0,28 


0,16 


0,22 


















c 


1,5 


0,4 




0,5 


0,1 


0,1 


0,1 


0,1 










d 


















0,302 


0,135 


0,04 


0,038 



The given above theory allows to build the dependency of thermodynamic character- 
istics on the chemical potential /i (charge carrier density n) in the wide interval of its 
values. 

However, within the study of doped MgB2 is of interest to examine /j, near the value 
of /iq ~ 0, 74:eV, which corresponds to the case of the pure MgB2. 

As it follows from the band calculations [4], /iq lies in the range of values, close to 
the complete filling of energy band a. For certainty we select the parameters Ci = C2 = 
0, (^ci — 0, 8eV, = 1) OeK. The interval where two bands overlap on the Fermi 

surface ^2 < A* < d gives the main contribution to superconductivity. It is convenient 
to represent the dependency of thermodynamic values from the value d = (fi — /xq) 7/^0 = 
(n — ho) /ho, which determines the relative change of charge carrier density as a result 
of substitutional impurity introduction into MgB2. At 5 > the density of electrons 
increases, and at 5 < - the density of holes. This approach allows to compare the 
obtained results with the experimental data (see Fig. 4 in the work [13]), because from 
these data is possible to determine the relative change of charge density S. 
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In Fig. 1 the dependency of from the parameter S is given. This dependency 
was obtained on the basis of developed above theory (curve 1). The same figure shows 
the points, which correspond to the experimental values of for diverse substitutions 
of atoms of Mg and B, taken from Fig. 4 of the experimental work [13], and also the 
experimental dependency (curve 2). 

From Fig. 1 follows, that the superconducting transition temperature Tc diminishes 
with the increasing of electron density (5 > 0) and remains constant at the introduction 
of holes into MgB2 (6 < 0). 

The decreasing of quantity Tc with the increasing of 6 is confirmed by experimental 
data for MgB2-xCx and M gQ^Q^CuQ^Q^B2-xCx- For these compounds the parameter 5 
is determined by the relations 5 — {x — 0, 05)/8, in accordance with the valence of the 
elements, which belong to the named compounds. 

In the case of Mgi_xLixB2, Mgi_xCuxB2, MgQ sLio^2B2-xCx for the quantities S we 
have respectively 5 = —x/8, —x/8, (—0, 2 + x)/8. This theory gives a constant value for 
critical temperature Tc — 39, AK with the condition of variable and negative 5 which takes 
place at any values of x in the case of two first compounds and at x < 0.2 in the case of 
the last compound. Figure 1 demonstrates the qualitative agreement between our theory 
and the experiment. This agreement is based on the idea of two bands and energy bands 
filling factor. The theoretical curve lies above the experimental points. This difference 
is related with two factors: firstly, the theory parameters were determined on the basis 
of the value Tc — 39, for the pure MgB2, (in the experimental work [13] Tc — ?>^K); 
secondly, the interband scattering of electrons on impurity was not considered. 

The theoretical dependencies of energy gaps Ai and A2 on 8 parameter in the domain 
of values which are of interest to compare with experimental data arc shown on Fig. 2. The 
values of these quantities for pure MgB2 correspond to Ai = 7, OSmeV, A2 = 2, Q8meV , 
remain constant when the system is doped with holes, and decrease with the increasing 
of electrons number in the system. Thus, the behavior of these quantities as a function 
of parameter is analogous to the behavior of the quantity T^. 

On the Fig. 3 the dependencies of the ratios ^i/KbTc (curve 1) and 2/S.2/ KbTc 
(curve 2) from relative change of charge density 5 are presented. For pure MgB2{S = 0) 
we have 2Ai/KbTc — 4.18, 2/^2/KbTc — 1.58. These values remain constant at 5 < 
and decrease at 5 > 0. 

Fig. 4 The relative jump of electronic specific heat (C5 — Cjv)/Cjv in the point T — Tc 
as a function of parameter 5 in MgB2- 

The relative jump of electronic specific heat in the point T = Tc as a result of of MgB2 
doping is shown in Fig. 4. As it follows from this figure, at 5 < the doping leads to 
a weak increasing of this quantity and to a much stronger increasing at 5 > 0. As was 
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expected, if the value of the jump (C5 — Cn)/Cn — 1.43 at T = Tc, that corresponds to 
the presence of one energy band. 

For certainty the dependency of chemical potential /i from the parameter 6 in MgB2 
is built on the basis of the second equation from (24) and is shown on the Fig. 5. This 
relation made possible to pass from the dependency of the given above thermodynamic 
quantities from chemical potential /i to charge carrier density, and also to the relative value 
S = .{n — nQ)/nQ Such transition allows comparing our two-band theory with experimen- 
tal data concerning the substitution of Mg and B by other elements of the Mendeleev's 
periodic table. 

5. The scattering of electrons on impurity potential. 



In the previous sections, on the basis of phonon mechanism, the theory of superconduc- 
tivity in two-band systems with reduced charge carrier density and nonmagnetic impurity 
was built. In this case the impurity influence is determined by the change of chemical 
potential /x (charge carrier density), which leads to the manifestation of the energy bands 
filling effect. The agreement with the experiment can be improved, if we will consider the 
processes of intra- and inter-band scattering of electrons on impurity potential. These 
processes affect not only the value of chemical potential /i, but also the thermodynamic 
characteristics of two-band system [9], [19], [21]. 

Further we will base on the results of the study [25] where the problem of the influence 
of nonmagnetic impurity on the superconducting transition in two-band systems with 
reduced charge carrier density is solved. With the precision of terms linear on impurity 
concentration we have: 



no + l/ni^Yl ^rn [Ccm + Xm - jJ^ - {Cm + Xm - Ij)- 



\Ccm ~l~ -^m A*| ~l~ I Cm ~l~ -^m A*| ~l~ Vm^^ i ^ 



(48) 



where Uq is the charge carrier density of pure substance, u - the difference between the 
valence of introduced atoms and atoms of the main substance. The quantities and 
Um are determined with the precision of terms linear on impurity concentration by the 
following relations: 



Xm. — 



^ ' ' n ' m.n. 



Cc 



E 



2r„ 



(49) 
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As it follows from the expression (48) the chemical potential /i depends explicitly on 
impurity concentration and through the parameters of which are determined by the 
relations (m, n = 1,2): 

^ = niTImn (50) 

where fii — ni/2Ni, and rjmn is determined by the inter- (n ^ m) and intraband {n — m) 
potentials of electrons scattering on impurity [25] which depend on Bloch functions and 
of unit cell volume. 

For superconducting transition temperature we obtain 

rc = Tco-«^^(l + ^) (51) 

0Ti2 iVi 

where Teg is determined by the formulas (28) and (29) taken in the respective value 
intervals of chemical potential (intervals 2 and 3) related with Ui through the relation 
(48), and a is determined by the expression [9], [19], [25]: 



Ni- N2 - - Ni - ^2 - 



X 



(All - A22)' - 4A12AJ "^\. (52) 



Here <1 (the sign + or — is selected considering the maximal value of Tc, which 
ensures a steadier superconducting state in comparison with the normal one). The ex- 
pression (51) is correct under the condition (2r„m7rTc)~^ <C 1, which takes place at low 
values of impurity concentration rii or at parameters r]njn- 

In Fig. 6 the dependencies of critical temperature Tc on the relative charge carrier 
density 5 obtained on the basis of expressions (51) are shown. Curve 1 in this figure 
corresponds to the absence of electron scattering on impurity {r]nm = 0,n;m = 1,2) and 
coincides with curve 1 in Fig. 1. Curve 2 corresponds to the values of 7711 = 0, 65, r/21 = 
0, 09. The dashed curve is an experimental dependency [13]. As it follows from this figure, 
the assumption of filling mechanisms for two energy bands and the electron scattering on 
impurity potential leads to a good agreement between theory and experiment in the 
case of systems doped with electrons (compounds MgQ_Q:^Cu(}Q Q^B2xCx and MgB2-xCx)- 
Substituting the atoms of Mg by Li and Cu we do not consider the electron scattering 
on impurity, assuming that in layered compounds of MgB2 the basic contribution to the 
superconductivity is given by the charges related to boron. The substitution of boron by 
C atoms leads to the violation of lattice periodicity in this layer and to the appearance 
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of scattering processes on impurity. However, the substitution of Mg by the atoms of 
Cu or Li (compound of Mgi_xLixB2, Mgi_xCuxB2) only changes the effective valence of 
boron, which is responsible for superconducting state. 

This idea is confirmed by the 5 dependency of superconducting transition temperature 
Tc on 5 in the compounds Mgo sLiQ^^B^-xCx. In Fig. 6 this dependency is characterized by 
the curve 3, which passes near the experimental points for this compound. The decreasing 
of Tc along with the decreasing of holes is obliged to charge scattering on impurity due to 
the substitution of boron by C atoms. The curve 3 in this figure is built on the basis of 
expression (51), in which the value of T^, equal to 39, iK is undertaken and the parameters 
of Vn — 0) 5; 7^21 = 0, 06 are selected. The difference between these parameters and the 
parameters inherent for MgB2-xCx, possibly, is caused by the difference between the unit 
cell volumes of these compounds [13]. 

Let's note, that during the construction of this theory the influence of doping on 
electronic state density A^i and A^2, and on characteristic phonon frequency luq wasn't 
considered. The consideration of this influence probably can change the values of the 
parameters of r)nm- 

6. Conclusions 



The theory of thermodynamic properties of two-band superconductors with variable 
charge carrier density on the basis of phonon mechanism of superconductivity was pre- 
sented in this work. This theory can describe the behavior of such values like Tc, Ai, 
A2, {Cs — Cn)/Cn at T = Tc as functions of charge carrier density in MgB2 when Mg 
and B are substituted by other elements from Mendeleev's periodical table. 

For this aim is necessary to do the following actions: 

1. To proceed from the BCS-type equation for the two-band model [7] with electron- 
phonon interaction constants, renormalized due to strong electron-phonon interaction and 
Coulomb interaction of electrons. 

2. To examine the interval where two-dimensional a and three dimensional tt bands 
overlap on the Fermi surface. To build the dependencies of the above-mentioned quantities 
from the chemical potential taking into account that the value = //q ~ 0, lAeV in 
MgB2 is near to the upper edge of a band, which is responsible for superconductivity [4] . 

3. To introduce the relative charge carrier density 5 = {fi — Ho)/l^o, which coincides 
with corresponding value calculated on the basis of the valences of elements, belonging to 
the compounds Mgi^xLixB2, Mgi^xCuxB2, Mgo,8Lio,2B2-xCx, Mgo^g5Cuo,o5B2-xCx and 
MgB2-xCx at different values of x. The S dependencies of Tc which are built in such a 



16 



way allow to compare theoretical and experimental data (see Fig. 1). 

The theoretical results obtained in this work qualitatively agree with the experiments 

results [13]: 

a) Namely, the introduction into MgB2 of electrons (5 > 0) decreases the critical 
temperature Tc{MgB2-xCx and Mg'o,95C'iio,o5-B2-xCx). 

b) The doping with holes {5 < 0) does not change the Tc {Mgi_^LixB2, Mgi_xCuxB2) 
with the change of 5. 

c) The compound MgQ sLio,2B2-xCx [13] [13], in which reaches the value of 39, AK 
when 6 = —0, 02 (that corresponds to MgB2) and diminishes with hole density decreasing, 
can't be framed in this scheme. 

5. The obtained results demonstrate also the quantitative agreement of the theory 
proposed in this work with experimental data for superconducting transition temperature 
Tc (see Fig. 6). 

a) In this case together with energy bands filling factor (or change of chemical po- 
tential /i) is necessary to consider the scattering of electron on impurity potential when 
the C atoms replace in the laminar structure the atoms of B, which is responsible for 
superconductivity. 

Upon consideration of these two mechanisms we obtain a dependency (curve 2), well 
describing the experimental data {MgB2-xCx and M5fo,95C''Uo,o5-B2-a;C'x)- 

b) Doping with holes 6 < {Mgi_xLixB2, Mgi_xCuxB2) does not change the value 
of Tc, since into the layer, responsible for superconductivity, impurity is not introduced, 
and the elements Li and Cu introduced instead of Mg only change the effective valence 
of S. 

c) The decreasing of Tc value (curve 3) in M go^sLio,2B2-xCx is related with electron 
scattering on impurity potential of C atoms. 
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